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We investigate the current-voltage (IV) characteristics of a model single-electron transistor where 
mechanical motion, subject to strong dissipation, of a small metallic grain in tunneling contact with 
two electrodes is possible. The system is studied both by using Monte Carlo simulations and by using 
an analytical approach. We show that electromechanical coupling results in a highly nonlinear IV- 
curve. For voltages above the Coulomb blockade threshold, two distinct regimes of charge transfer 
occur: At low voltages the system behaves as a static, asymmetric double junction and tunneling 
is the dominating charge transfer mechanism. At higher voltages an abrupt transition to a new 
"shuttle" regime appears, where the grain performs an oscillatory motion back and forth between 
the electrodes. In this regime the current is mainly mediated by charges that are carried on the 
grain as it moves from one electrode to the other. 



I. INTRODUCTION 

The mechanical properties of mesoscopic conductors 
and their influence on charge transport are very much in 
the focus of recent solid state physics research. Certain 
anomalous behavior of nanowireada and the electrostat- 
ically controlled deformation of carbon nanotubesci are 
examples of an interplay between electrical and mechan- 
ical degrees of freedom that appear on the nanometer 
length scale. Other examples where the heteroelastic na- 
ture of a material crucially affects single-electron tunnel- 
ing have been found in studies of self-assembled metal- 
organic composite structures. 

The relevant scenario associated with a strong elec- 
tromechanical coupling is that significant deformations 
occur as a result of large Coulomb forces acting on charge 
accumulated in some small region, for instance, in a 
metallic cluster. Recently, a model system containing 
such a coupling was considered by Gorelik et alu, who 
proposed a single-electron tunneling device containing a 
movable metallic cluster in tunneling contact with bulk 
metallic electrodes. In this device mechanically soft or- 
ganic links serve both as elastic springs, keeping the 
cluster in place, and as tunnel barriers with resistances 
that are exponentially sensitive to the deformation of the 
springs. An important consequence of the interplay be- 
tween single electron tunneling and the mechanical vi- 
bration of the cluster in this model is the electromechan- 
ical instability predicted in Ref. |^: If a large enough bias 
voltage is applied between the electrodes, the equilibrium 
position of the grain looses its stability and cluster vibra- 
tions develop. Such vibrations give rise to a new mecha- 
nism of charge transfer, where electrons are transported 
through the system by the metallic cluster which per- 
forms shuttle motion between the electrodes. The elec- 
tric current, / = 2Nef , associated with this mechanism 
does not depend on the tunnel transparencies, and is only 
determined by the frequency, / = uj/2n, of the elastic vi- 
brations of the cluster and the number, N, of electrons 



carried by the cluster. Experimental evidence for a cou- 
pling between electron transfer and vibrational degrees 
of freedom iias been found both for macroscopidj and 
microscopicOij systems. Different aspects of this phe- 
nomenon has, also been theoretically investigated in sev- 
eral articlesEJlla |-||-| 

In the work discussed aboveou it was shown that a 
large damping constant, 7, is detrimental for the devel- 
opment of the shuttle instability and in the limit where 
7 > /, elastic shuttling of the charge becomes impossi- 
ble. The mechanical lability of the system, however, is 
still a dominating feature of the charge transport even in 
the limit of strong dissipation. What the consequences 
of such a lability would be is a question which needs to 
be answered. This is not only an academic question since 
coupling to intramolecular vibrations inside deformable 
organic molecules carrying current as well as friction in 
the medium through which a metallic cluster is moving 
may cause significant dissipation. The dissipative limit 
of electromechanical mesoscopic structures with movable 
parts is therefore important for understanding the func- 
tioning of realistic nanometer size structures. Our objec- 
tive in the present work is to study this limit. 

We will consider the model system illustrated in Fig. |^. 
The current flow between the metallic electrodes is due 
to the tunneling of electrons between the electrodes and 
the metallic cluster. This is assisted by the displacement 
of the cluster. An electrostatic force acts on the charged 
grain if a finite bias voltage is applied between the elec- 
trodes. The one-dimensional dynamics of the cluster is 
also governed by an elastic restoring force and a friction 
force. In contrast to the approach developed in Ref. ^ we 
will consider the limit where the electric force dominates 
over the elastic force, which means that the dynamics 
of the charged cluster is determined by the interplay be- 
tween Coulomb forces and friction. This, however, does 
not mean that the elastic forces can be totally neglected. 
For bias voltages slightly above the Coulomb blockade 
threshold when the cluster is in the center of the system. 
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the dynamics of the cluster is actually very sensitive to 
the value of the elastic force. The low temperature non 
linear charge transport through the system is affected 
both by the Coulomb blockade phenomenon and the me- 
chanical motion of the cluster. These two phenomena are 
coupled since the threshold voltage for electron tunnel- 
ing depends on the junction capacitances which in turn 
depend on the cluster position with respect to the elec- 
trodes. A general property is that the threshold voltage 
increases when the distance between the cluster and an 
electrode decreases. 

In order to understand qualitatively the electrome- 
chanical charge-transport scenario, let us consider a neu- 
tral cluster situated in its equilibrium position between 
the electrodes where the voltage threshold for electron 
tunneling has a minimum value, Vq. At zero tempera- 
ture no tunneling is possible for voltages V < Vq, where 

V is the bias voltage applied between the electrodes. For 

V > Vo tunneling onto the cluster becomes possible and 
the cluster can be charged. It is easy to understand that 
the direction of motion of the charged cluster, due to the 
Coulomb forces, will be away from the electrode which 
has supplied the extra charge to the cluster. After some 
time the extra charge will disappear, usually to the near- 
est electrode, which makes the cluster charge zero again. 
An important question at this stage is whether one more 
tunneling event to the nearest electrode is possible or not. 
The answer is not evident since the electrostatic thresh- 
old is different from the one at the initial point in the 
middle of the system. As we will see, depending on the 
applied bias voltage, we can have one of two possible situ- 
ations. For voltages Vq < V < Vt, where Vt is a thr eshold 
voltage which will be treated in more detail in Sec. IV C , 



the extra tunneling event is not possible. In this case the 
cluster is almost trapped near the electrode. Small oscil- 
lations in the vicinity of the trapping point are possible 
due to the action of the weak elastic force, but the clus- 
ter will not be pushed back by Coulomb forces. In this 
case the conductance is not assisted by significant clus- 
ter displacements between the electrodes. We call this 
regime the tunneling regime since the charge transfer is 
very similar to the conventional single electron transport 
in a static system. 

li V > Vt there is a possibility for another tunneling 
event between the grain and the nearest lead to happen 
after the extra charge has tunneled off the cluster. This 
event changes the sign of the net charge on the grain. 
In this case the cluster can be pushed by the Coulomb 
force towards the more distant electrode where the above 
described process repeats itself. The conductance is now 
assisted by significant displacements of the grain and this 
scenario is qualitatively similar to the shuttle vibrations 
in fully elastic electromechanical structureafl. We call this 
regime the shuttle regime of charge transport. A sharp 
transition, corresponding to a current jump, occurs in a 
small voltage interval between the two regimes. 



II. MODEL SYSTEM 

We will consider a model based on the picture in Fig. [|. 
This is a simplified model which, however, retains many 
of the interesting features of a "real" system. The system 
consists of a metallic grain of mass M placed in the gap 
between two bulk leads separated by a distance L. The 
displacement of the grain from the center of the system 
is measured by the coordinate We consider only ID- 
motion of the grain between the leads. A bias voltage V 
is applied between the leads. In this simplified case we 
take into account only three different forces acting on the 
grain; a linear elastic restoring force F^i = —kX, a dissi- 
pative damping force Fdiss — —IdX, and an electrostatic 
force Fq. The electrostatic force is a function of the bias 
voltage, V, and the charge, Q, on the grain: 
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Here Cq is a capacitance constant determined by the ge- 
ometry of the system. To get this expression we assume 
that the capacitances Cl between the left lead and the 
grain and Cr between the grain and the right lead can 
be approximated as parallel plate capacitors and that all 
other capacitances can be neglected. The first term in 
Eq. (Q) can be understood as the force from an effective 
electrostatic field V/L in the junction, which couples to 
the extra charge on the grain. The second term can be 
thought of as the interaction of the charge on the grain 
with image charges in the two leads. Note that the last 
term in Eq. (0) always results in an attraction of the 
charged grain towards the nearest lead. If we take these 
forces into account we can write the equation of motion 
for the grain as: 



QV 

MX + jdX + kX = — 



X 



(2) 



We can now consider the Coulomb blockade regime 
where the Coulomb charging energy, — e^/2C is larger 
than both quantum and thermal fluctuations, 3> 
h/ RC, /3~^ , where R is the smallest tunneling resistance 
possible in the system, and /? is the inverse temperature. 
We thus assume that 



RiX) » Rq 



2? 



6.5kn. 



for all positions X available for the grain^. If the cri- 
teria for the Coulomb blockade regime are met, we can 



*Our approach is based on a classical description of the 
grain displacement and is different from approaches where 
quantum akister vibration assisted tunneling to the grain is 
consideredtJ 

^Note that soft matter springs always have some finite thick- 
ness even when compressed. Coating layers on the leads could 
also restrict the space available for the grain. 
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consider electrons on the grain to be fully localized and 
express the charge on the grain as 

Q{t) = en{t), 

where n{t) takes on only integer values, (e is the electron 
charge.) Let (n, Qa) be the state of the system with n ex- 
tra charges on the grain and the charge Qa on the lead a 
(a = L, R). It then jpUows from the "orthodox" Coulomb 
blockade theoryll^O that the tunneling probabilities for 
the tunneling event (n, Ql,r) {n^l, Ql,rTs) to occur 
during the time At are 

Pij,in,X,V,At) = 

AGl^in,V,X) r _ / AGt,^in.V.X) \-] ~' 

(3) 

where AG'^^(n, X) is the decrease of free energy in 
the system as an electron tunnels to the right (+) or 
to the left (— ) through the left {L) or right [R) tunnel 
junction, ks is the Boltzmann constant, and Rl,{r.) is the 
resistance of the left (right) tunnel junction. This resis- 
tance depends exponentially on the displacement of the 
grain from the center of the system and can be written 

RLiX) - Rr{-X) = i?oexp(X/A), 

where Rq is a constant prefactor and A is referred to as 
the tunneling length. The tunneling length depends on 
the materials used in the system and for our system we 
estimate A to be of the order of one Angstrom. 

III. CURRENT- VOLTAGE CHARACTERISTICS 

To make the treatment of the model easier we rewrite 
the equations in a dimensionless form. If we introduce 
the dimensionless time r = t/tg where to — Id.L'^ /gVq 
is a timescale on which the grain crosses the distance 
L between the leads due to the electrostatic forces, the 
dimensionless length x — X/ the dimensionless elastic 
vibration frequency lo = -^//cL^/eVo, the dimensionless 
bias voltage v = V/Vq where Vg = e/4Co is the Coulomb 
blockade threshold in the center of the system, and the 
dimensionless constant a — MeVb/i^7^ which signifies 
the ratio between the electrostatic force and a typical 
dissipative damping force in the system, we can rewrite 
the equation of motion for the grain as: 

ax + X + Lip'x = nv + An'^x (4) 

We will focus on the case when the dissipative force dom- 
inates the electrostatic force while the latter dominates 
the elastic forces, ^ a ^ 1. The free energy terms 
AG^ to be used in Eq. (^) are 



AGt{n,v,x) = ^{l~ix')(^-lT2n±^^^ (5) 

AG±(n, v,x) = ^{l^ ix') (^-l±2n± . (6) 

Note that the position dependence of AG^ ^ given by 
Eqs. (H) and results in a position dependent Coulomb 
blockade threshold voltage. This means that whether 
tunneling in a junction is blocked or not at a certain 
voltage depends on where the grain is located at the mo- 
ment. 



A. Numerical approach 

In a numerical approach we have peformed Monte 
Carlo simulations of the model system described. A 4:th 
order Runge Kutta method was used to solve the equa- 
tion of motion for the grain for small enough time steps 
for the charge on the grain to be considered constant 
during each step. After each time step the charge on the 
grain was updated by "rolling dice" and deciding whether 
to carry out a tunneling event using the tunneling prob- 
abilities of Eq. (||). The current was calculated as the 
average of the number of transfered electrons over a cer- 
tain time interval. For our choice of parameters the av- 
erage number of electrons transfered through the system 
stabilizes over a time period of approximately 64^0- In 
our calculations of the current - voltage characteristics 
we have averaged over 6.4 • lO'^io = 20 /iS to reduce the 
numerical noise. The result of the calculation is plotted 
in Fig. |. 

If we now compare the current through the studied sys- 
tem with that through a static symmetric double junc- 
tion as in Fig. || (see inset), it becomes very clear that 
there are two distinct parts of the current-voltage curve. 
For voltages V approximately between Vq and 1.5Vo the 
current through the static symmetric double junction is 
larger than that through the system under considera- 
tion here. Since, as will be shown in the next section, 
charge transport in this regime is dominated by tunnel- 
ing, we label this regime the tunneling regime. For volt- 
ages above approximately 1.5Vb, on the other hand, the 
current through the present system is the larger one, and 
since, as will be shown below, charge transport in this 
regime is mechanically mediated by the grain, we label 
this regime the shuttle regime. 

The distinction between the two regimes also becomes 
very clear if we consider the root-mean-square of the dis- 
placement of the grain from the center of the system, x, 
as a function of the bias voltage. This is plotted in Fig. ||. 
It is clear that the average displacement is much larger for 
the tunneling regime than for the shuttle regime. The av- 
erage displacement is also increasing with the bias voltage 
for the tunneling regime, whereas, for the shuttle regime 
it is a slowly varying function. 
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IV. DISCUSSION 

A. Tunneling regime 

From Fig. ^ (inset) we see that for bias voltages just 
above the Coulomb blockade threshold Vq (for the grain 
in the center position) the current is smaller than it is 
through a static symmetric double junction. To under- 
stand this we should consider the x-dependence of the 
tunneling rates. Assume that the grain starts out sit- 
ting uncharged in the center of the system and that the 
bias voltage is just above Vq. At this point two things 
are possible. One unit of charge can either tunnel onto 
the grain or off the grain. Since the system is symmetric 
we consider only the first of these cases. The criterion 
for tunneling from the left lead to the grain is that the 
free energy is lowered after a tunneling event, AG J > 0. 
Using Eq. (^) we find the corresponding inequality 



Note that x is the normalized coordinate so that 
— 1/2 < X < 1/2. We see that the Coulomb blockade 
threshold when the grain is at the center of the system 
is (x = 0) is indeed Vq {v = 1). For lower voltages tun- 
neling onto the grain from the left lead is still possible 
as long as the grain is to the right of the center position. 
However, this process is exponentially suppressed due to 
the increase of resistance with tunneling distance. If one 
considers tunneling from a neutral grain to the right lead 
the same picture (with x —x) emerges. When the bias 
voltage is increased above Vq, tunneling onto the grain 
becomes possible if it is to the left of the center. We see 
here that if the bias voltage is not much higher than the 
Coulomb blockade threshold Vq, the open region, where 
both tunneling onto a neutral grain from the left lead 
or off a neutral grain to the right lead is allowed at the 
same time, is much smaller than the distance between 
the leads. The concept of the open region is illustrated 
with two examples in Fig. |4[ 

For a grain that has the charge n = and is located 
inside the open region, both the processes n +1 and 
n ^ — 1 are allowed at the same time. If the grain is 
located outside the open region it can only be charged 
from the far lead. 

Let us now consider the case when the bias voltage is 
not much higher than the Coulomb blockade threshold Vq 
that applies if the grain is in the center position. In this 
case the open region is much smaller than the distance 
between the leads. If a unit charge tunnels onto the grain 
from the left, the grain becomes positively charged and 
is thus affected by a force towards the negative (right) 
lead. It will start to accelerate towards that lead, but 
if the mass of the grain is very small and the dissipa- 
tion large, the grain will reach a maximum velocity very 
quickly. As the grain comes close to the negative lead. 



the decharging process through the right junction be- 
comes very probable. If the relaxation of the charge on 
the grain to the negative lead takes place outside the 
open region, the grain cannot be recharged by a nega- 
tive unit charge from the negative lead. If dissipation is 
strong the grain will stop very quickly and the very small 
elastic restoring force will start to move the grain very 
slowly towards the center of the system. At this time the 
grain is only in tunneling contact with the far lead and 
it will continue to move slowly towards the center, either 
until it reaches the open region and can be charged from 
either lead or until a tunneling event from the positive 
lead on the far side of the system occurs again. If the last 
of these two processes occurs, the charge on the grain be- 
comes positive and the grain is accelerated towards the 
negative lead again, repeating the above described pro- 
cess. The resulting motion is thus an oscillation around 
an average position, which is located between the open 
region border and the lead. One cycle of such an oscilla- 
tion is schematically illustrated in Fig. |[ 

Tunneling from the far lead to the grain as the grain 
moves under the influence of the weak restoring force is 
possible but very unlikely, as can be seen from Eq. (|^). 
This is due to the exponential dependence of the tunnel- 
ing resistance on the separation between grain and lead. 
If the grain moves very slowly, however, there may be 
enough time for the grain to be charged from the far 
lead before it reaches the open region. As the bias volt- 
age is increased, the size of the open region increases, 
thereby affecting the probability that the grain will reach 
the open region before getting charged from the far lead. 
This leads to a transition to the shuttle regime discussed 
in the next section. 

We can thus conclude that the current through the 
system is smaller than that through a static symmetric 
double junction because the charge transfer mechanism 
is limited by tunneling through the more resistive tunnel 
barrier, just as is the case for a static asymmetric double 
junction. That this is actually the case is also confirmed 
by studying plots of the grain position as a function of 
time, obtained from Monte Carlo simulations of the sys- 
tem. Such a plot for V — 1.1 Vb is shown in Fig. ||. For 
clarity the picture is embedded in a model system with 
the positions of the leads marked on the x-axis. The plot- 
ted line traces out the position of the grain as a function 
of time. The sharpness in the curve depends on the two 
very different time scales in the system. The time scale 
for grain motion due to the electrostatic force when the 
grain is charged is much smaller than the time scale of 
grain motion caused by the weak elastic force when the 
grain is uncharged. From the plot we can conclude that, 
on average, the charge transfer through the system looks 
like that through a static asymmetric double junction. 

If we consider Eq. (^ and its counterpart, 

AG+(n = 0,x) > 0, 

we see that for = l.lVb, the borders of the open region 
are located at a; = ±0.05. When we compare this value 
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to the average displacement of the grain at this voltage, it 
becomes clear that the average displacement is 3 - 4 times 
bigger. The grain thus performs an oscillatory motion 
around an average displacement, and these oscillations 
are possible because the average displacement is located 
quite far away from the border of the open region. 

We can now compare the current for this regime, (see 
Fig. |2|) with the current through a static asymmetric 
double junction. The current through the latter type of 
double junction can be approximated by saying that the 
charge transfer to the far lead limits the current. Since 
the inverse of the tunneling rate is the average time be- 
tween tunneling events, we can write the current la.d.j. 
through the asymmetric double junction as 



la.d.j. — f ar—te 



(8) 



where Tfar-iead is the rate for tunneling events between 
the grain and the far lead. Using Eq. under, for 

instance, the assumptions a; > and T = 0, the current 
from the far lead would be 



la.d.j. — 



SRqCo 



exp I 



(9) 



If we use the average displacement from Fig. I the cur- 
rent, as calculated by Eq. and in the voltage interval 
1 < V/Vq < 1.25, turns out to be of the order 20 % lower 
than the actual current through the our system. This is 
understandable since the small grain oscillations around 
the average displacement decrease the effective tunneling 
resistances seen by the charges transfered through the 
system. 



B. Shuttle regime 

The statements made in the previous section mean that 
we can expect the current through our device to be very 
small on the scale of the current through a symmetric 
static double junction. On this scale, we can also expect 
that the current only increases slowly as the bias voltage 
is raised to slightly above the Coulomb blockade thresh- 
old in the center of the system. The current will continue 
to increase very slowly with the voltage. As the size of the 
open region increases it becomes more and more proba- 
ble that an empty grain will reach the open region before 
it is recharged from the far lead. If the grain reaches the 
open region, charge transfer from the near lead suddenly 
becomes the dominating charge transfer mechanism. If 
we consider Eq. (0) we see that as the bias voltage V 
reaches 2Vb, the open region has extended all the way 
to the leads. The grain will thus always move inside the 
open region. In this case, the charge transfer cycle looks 
quite different from the picture in the previous section. 
When the grain gets positively charged it will move to- 
wards the negative lead. As the grain gets closer to the 
lead, the tunnel resistance decreases exponentially and 



finally the charge on the grain will tunnel from the grain 
to the lead. When the grain loses its charge it will stop 
very quickly due to the high dissipation. The grain now 
starts to move very slowly towards the center of the sys- 
tem, but since the timescale of charge exchange with the 
near lead is much shorter than that of movement due to 
the elastic force, another tunnel event can occur and the 
grain can get negatively charged. This means that the 
grain will be accelerated towards the positive lead, where 
a similar procedure will occur. The grain will now con- 
tinue to move back and forth in this fashion, shuttling 
charge across the junction. Figure ^ shows a schematic 
illustration of this charge transfer mechanism. 

The exponentially large tunnel resistances limiting the 
current in the tunneling regime are now gone, since all 
tunneling events occur when the grain is close to the 
leads. The oscillations of the grain thus effectively lower 
the tunnel barriers seen by the transfered charges, which 
leads to a large increase in the current. 

We can now proceed as in the case of the tunneling 
regime and consider a plot of the grain position as a 
function of time for some bias voltage in this interval. 
In Fig. ^ we have made such a plot for the bias volt- 
age V = 2.0Vb. As in Fig. ^ the plot is shown together 
with the model system so that the positions of the leads 
are marked on the x-axis. The plotted line traces out the 
position of the grain as a function of time. The grain per- 
forms a stochastic but still oscillatory motion back and 
forth through the system. For this voltage an uncharged 
grain is everywhere in tunneling contact with both leads 
so that the charge on the grain can change by 2e at each 
approach of a lead. This means that the grain will al- 
ways be pushed by the electrostatic force, which explains 
the much shorter time scale for grain motion in Fig. ^ 
compared to Fig. ^ (Note also the factor of 10 difference 
in scale on the time axes in the two plots). 

Let us now go back and consider the IV-cuTve in Fig. ^ 
again. For a bias voltage of approximately 4.5Vo, the 
IV-cvLTve changes slope over a relatively short voltage 
interval. The reason for this is the transition to a regime 
where two extra charges are allowed on the grain, i.e. 
four charges can be transported across the system in each 
shuttle cycle. To get a better understanding of this, we 
should consider the case AGj(n = l,a;) > 0, i.e. the 
condition that the free energy decrease should be posi- 
tive when one charge tunnels from the left lead onto an 
already charged grain. Using Eq. (^), we get the condi- 
tion 



X > 



1 



(10) 



We thus see that when V — 3Vb {v — 3) a. new open 
region develops, where electron tunneling is allowed from 
the left lead when the grain charge is n = 1 and to the 
right lead when the grain charge is n = — 1. If we remind 
ourselves of what went on in the tunneling regime, we 
cannot expect that the current will change much until 
the size of this region is of the same size as the amplitude 
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of the grain oscillations. As can be seen from Fig. ||, 
nothing new happens to the /^-characteristics when V = 
SVq. However, approximately when V = 4.5]^ there is 
a transition to the new regime. From Eq. ( |lO|) we get 
that, at V — 4.5, the open region borders for n — 1 have 
extended to approximately x € (—0.25,0.25). At this 
voltage, the grain oscillations should thus be inside the 
new region most of the time, allowing the transfer of four 
charges in each shuttle cycle. 

It is important to note here that the transition in the 
IV-cuTve is not sharp. As the open region for n = 1 grows 
wider, it will become more and more probable that, as 
the grain moves across the system, it will transport two 
charges instead of only one charge. As is normally the 
case for shuttle transporta, we can consider a current 
frequency relationship: 

/ = 2AAe/, 

where Af is defined by this equation and represents an av- 
erage number of extra electrons transported on the grain, 
and where / is the vibrational frequency of the grain. 
Both J\f and / are functions of the bias voltage. Note 
also that Af is not usually an integer. 



C. Analytical description of the shuttle regime 

In this section we present an analytical approach to 
modeling the current through the system for bias volt- 
ages in the range 2 < V/Vq < 3. In this voltage inter- 
val, the grain can only shuttle one charge at a time in 
each direction. Since the motion of the grain is strongly 
influenced by the random tunneling events, we have to 
consider the period time in an averaged sense and write 
the current as: 



toT 



(11) 



where T is the dimensionless average oscillation period 
and to = JdL'^ /gVq is the typical time scale in the system. 
We make the assumption that we can divide the average 
period into the three parts schematically illustrated in 
Fig. |9[ Since the system is symmetric with respect to 
the center of the system it is enough to consider half a 
cycle. 

The first part of the average period, Ti, is the aver- 
age time it takes a grain with one excess charge to move 
from the center of the system towards the negative lead 
to the position where, on average, the excess charge is re- 
laxed to the negative lead. After the charge has relaxed 
to the negative lead, the grain stops very quickly and, on 
the average, sits still during the time T2 before one more 
charge tunnels to the negative lead. As this happens, it 
takes the grain the time T3 to get back to the center of 
the system, where it repeats a mirror version of this cycle 
towards the positive lead. Note also that the further the 



grain moves towards the lead, the shorter the time T2 
can be expected to be. As the grain passes the center of 
the system towards one lead, there is at each position a 
certain probability that the charge on the grain will tun- 
nel to the lead. Wherever the tunneling event occurs, the 
two average times T2 and T3 are determined by the first 
time Ti, which is determined by the position at which 
the tunneling event occurred. We can therefore write the 
average period time as: 



T 



T{x)P{x)dx, 



(12) 



where t(x) is the half-period for a grain that reaches 
position X as it travels from the center of the system 
towards the lead. This half-period now consists of the 
sum of three partial times ti{x), T2{x) and T3{x), where 
the indexes refer to the same parts of the half-period as 
the time indexes illustrated in Fig. ^. 

To find the probability density P{x), we can consider 
an ensemble consisting of N grains. These grains all start 
out at the center of the system, have charge n — 1 and 
move towards the negative lead. We can first find the 
relative number of grains m{x)/N that still has a charge 
of n = 1 at X by noting that: 



dx 



m{x) rjj(n = 1, x) 



(13) 



This is an ordinary separable differential equation with 
the solution: 



m{x) m(0) 



N 



N 



■ exp 



- tor+{n^l,x'] 
x{x') 



dx' 



(14) 



Since all grains in the ensemble have charge n = 1 at 
T = 0, we see that ■m{0)/N = 1. We can now find the 
probability density P{x) as the relative number of grains 
in the ensemble that stops at precisely x, i.e. minus the 
derivative of m{x)/N; 



P{x) 



tpT^jn = l,x) 
x{x) 



exp 



l,x') 



x{x') 



dx' 



(15) 



The next step is to find the half-period, t(x). Since we 
are working in the high dissipation limit, a ^ 1, acceler- 
ation times are very short compared to the time scales of 
movement of the grain and tunneling. This means that, 
we can to a good approximation find the parts ti(x) and 
T3(a;)|-by integrating the equation for the velocity of the 
grairO: 



i = nw + 4n x, 



(16) 



from T = to r = t'{x), and for n = ±1. The resulting 
traveling times are: 



6 



Tl{x) 



■ In 1 



(17) 
(IS) 



To find the time T2 we first assume that the grain wiU 
not move on the scale of the tunnehng length during this 
time. This means that the tunneling rates are time inde- 
pendent and that we, if the grain sits with zero charge at 
X, can expect the average time before a tunneling event 
occurs to be: 



T2(x) 



1 



At zero temperature we can expect the time to be: 
8i?oCo exp(-fx) 



T2{x) 



to 



(l-4x2)(-l + TTk) 



(19) 



(20) 



We have thus arrived at the following expression for the 
current through the system in the bias voltage interval 
2 < V/Vn < 3: 



to Jn"'°" (Tlix) + T2{x) + T3(x)) P{x)dx' 



(21) 



where ti(x), T2{x) and t^^x) are given by the equations 
(0), (|6|) and (^) and P{x) is given by Eq. (|l|). 

We have, with the same parameters as used in our 
earlier Monte Carlo simulations, numerically calculated 
the current given by Eq. (21). The results are shown in 
Fig. 
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The solid line corresponds to the Monte Carlo 
simulations of the system and the circles correspond to 
the values obtained from Eq. (|l]). The agreement be- 
tween the numerical studies and the analytical approach 
is very good, which is a strong indication that the charge 
shuttle mechanism description of the charge transfer is 
applicable also in highly dissipative systems. 

It is also of interest to know the threshold voltage, Vt, 
and the width, AV^, of the transition from the tunnel- 
ing regime to the shuttle regime. In order to estimate 
these we can consider small oscillations. Ax, of the grain 
around some average position xq ■ Without loss of gener- 
ality we can assume that xq > 0, i.e. the grain oscillates 
on the right hand side of the system. If we assume that 
the oscillation amplitudes are not very big we can esti- 
mate the velocity of the grain to be 



VIN 



-uj'^xo, 



(22) 



for grains moving towards the center of the system due 
to the elastic force. On average it moves during the time 



1 



Tin 



to^ti'^ = 0,v,xo) 



(23) 



before it is charged from the far lead. When the grain 
is moving towards the lead due to the electrostatic force 



acting on the extra charge on the grain, it approximately 
moves with the velocity 



VOUT « (w + 4xo). 



(24) 



The average time it will move before the extra charge 
tunnels to the right lead is 



Tout 



toT]^{n = l,v,xo) 



(25) 



For the position xo to be stable the average the distance 
the grain moves in each direction has to be equal to each 
other. We thus get the relation: 



OJ Xq l 

-e>- 



h 



.0 _ t±^e-T-o 



(26) 



where /l = -1/2 + CrV/c and = 1/2 Ci^/e are 
functions of the right and left capacitances and the bias 
voltage and that are of order unity as long as the grain 
is not close to the open region border. Rearranging the 
factors in Eq. (^6|) and taking the logarithm of both sides 
we get 



^^Xq = In^ 



ln!i±i^+ln4«. 



Xq 



(27) 



Under the conditions that we are not close to the open 
region border and that the elastic force is very weak we 
can neglect the last two terms on the right hand side of 
Eq. ( P7| ) . In this case we get the average position for the 
grain as 



Xo 



A , 2 



(28) 



If Aa; ^ Xq, one can, by comparing the average position, 
Xo, for the grain with the open region border, (v — l)/2, 
estimate the threshold voltage. 



(29) 



which corresponds to the transition from the tunneling 
regime to the shuttle regime. 

We can now use the expression for xq to estimate the 
width of the oscillations as 



Ax 



X RoCo r^. 1 
2 In 

L to iu 



(30) 



From Eq. (|^) we know that the open region expands 
linearly with the bias voltage. When the oscillations 
are completely outside the open region we can expect 
the system to be in the tunneling regime. When the 
open region has expanded to include the oscillations, the 
system should be in the shuttle regime. The open re- 
gion border expands Ax if the voltage is increased with 
AV/Vo = 2Ax and we thus get the relative transition 
width as: 



7 



Vt-Vo 



Ax 



4 77 



(31) 



where Ws/;. — ^o^'^ = cVo/^dL'^ is a typical grain oscillation 
frequency, loj^ = 1/ RqCq is a characteristic tunneling fre- 
quency, and rj = Lo"^ ~ kL'^/eVo represents the strength 
of the electromechanical coupling. From Eq. (^ij) one 
can see that there are two cases when there is a very 
sharp transition between the two regimes. The first case 
is when the electromechanical coupling becomes very 
strong. The second case is when the shuttle frequency 
is low compared to the rate of tunneling. In our system 
these conditions arc realized by the assumed weak elastic 
forces and the high rate of dissipation associated with the 
moving grain. 



V. CONCLUSIONS 

The main conclusion resulting from our analysis is that 
an electromechanical coupling in dissipative nanometer 
sized Coulomb blockade structures cannot be viewed sim- 
ply as an additional channel for absorbing the power as- 
sociated with the current injected into the system. In- 
stead a new mechanism of mechanically assisted charge 
transfer occurs, which increases the current exponentially 
and which to some extent is related to the shuttling of 
electrical charges, predicted for weakly dissipative elec- 
tromechanical structuresQu. We have shown that the 
electromechanical coupling results in a highly nonlinear 
IV-curve with two distinct regimes of charge transport. 
More features of the charge transfer might be available 
by studying the noise properties of the system. Since the 
noise is sensitive to the dynamical properties of the sys- 
tem, noise measurements can give additional information 
about the interplay between elasticity and dissipation in 
real nanoelectromechanical structures. 



* L.Y. Gorelik, A. Isacsson, M.V. Voinova, B. Kasemo, R.I. 

Shekhter, and M. Jonson, Phys. Rev. Lett. 80, 4526 (1998). 

A. Isacsson, L.Y. Gorelik, M.V. Voinova, B. Kasemo, R.I. 

Shekhter, and M. Jonson, Physica B 255, 150 (1998). 

M.T. Tuominen, R.V. Krotkov, and M.L. Breuer, Phys. 

Rev. Lett. 83, 3025 (1999). 
^ A. Erbe, R.H. Blick, A. Kriele, and J.P. Kotthaus, Appl. 

Phys. Lett. 73, 3751 (1998). 
® A. Erbe, C. Weiss, W. Zwerger, and R.H. Blick, Phys. Rev. 

Lett. 87, 096106 (2001). 
^H. Park, J. Park, A.K.L Lim, E.H. Anderson, A.P. 

Alivisatos, and P.L. McEuen, Nature (London) 407, 57 

(2000). 

^° C. Weiss, W. Zwerger, Europhys. Lett. 47, 97 (1999). 
" D. Boese, H. Schoeller, Europhys. Lett. 54, 668 (2001). 

N. Nishiguchi, Phys. Rev. B 65, 035403 (2001). 

L.Y. Gorelik, A. Isacsson, Y.M. Galperin, R.I. Shekhter, 

and M. Jonson, Nature (London) 411, 454 (2001). 
" A. Isacsson, Phys. Rev. B 64, 035326 (2001). 

D. Fedorets, L.Y. Gorelik, R.I. Shekhter, and M. Jonson, 

cond-niat/0104200| (unpublished). 

R.I. Shekhter, Sov. Phys. JETP 36, 747 (1973). 

I.O. Kulik and R.I. Shekhter, Sov. Phys. JETP 41, 308 

(1975). 

Note that P^^ = ^t^L,R- 
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FIG. 1. Schematic picture of the model system, which consists of a mctalhc grain of mass AI coupled by weak clastic links 
to two electrodes separated by a distance L. The elastic links act as springs with spring constant k. The tunneling resistances 
of the right and left junctions are Rr and Rl- A bias voltage V is applied across the system. 



9 




1 2 3 4 5 6 



V/Vo 

FIG. 2. The solid line shows the current - voltage characteristics obtained by a Monte Carlo simulation of the charge transport 
through the system sketched in Fig.^. The calculated current, which was averaged over 20/is, is plotted as a function of the 
bias voltage V scaled by the Coulomb blockade theshold voltage Vq, that applies if the movable grain is equally far from both 
electrodes. The dashed line displays the current through a static symmetric double junction for the same parameters. The 
parameters used in the simulation are: a — 6.4* 10~* and uj^ = 4.27* 10"'^. It is clear that for voltages between approximately 
Vo and l.SVb (see the inset which shows a magnification of the voltage interval 1 < V/Vo < 2) the current through the model 
system is smaller than the current through the static symmetric double junction, whereas, for higher voltages it is the other 
way around. 



10 




°0 1 2 3 4 5 6 

V/Vo 

FIG. 3. The root-mean-square displacement of the grain from the symmetric position between the leads as a function of the 
bias voltage scaled by Vb, the Coulomb blockade threshold voltage in the center of the system. The parameters used in the 
simulation are: a = 6.4* 10"'' and ijp' = 4.27* 10~^. The distinction between the two different regimes of charge transfer is very 
clear. In the tunneling regime, the average displacement increases with the voltage and is larger than in the shuttle regime. In 
the shuttle regime, the average displacement is a slowly varying function of the voltage. 
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(b) 




FIG. 4. Illustration of the concept of the open region which, in the pictures above, correspond to the space between the 
vertical solid lines, (a) If the grain is uncharged and located inside the open region, it is in tunneling contact with both leads at 
the same time, (b) When the grain is situated outside the open region, energy considerations show that tunneling to the near 
lead is blocked. Tunneling from the far lead is still possible, however, this process is strongly suppressed due to the exponential 
dependence of the tunneling resistance on the grain-lead separation. 
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1: The grain is charged from the far lead. 




I 



2: Being charged, the grain is pushed towards the 
near lead until a tunnel event to that lead occurs. 



+1 




3: The grain is slowly pulled back by the weak elastic 
force until charging from the far lead occurs again. 

FIG. 5. Schematic illustration of the charge transfer mechanism in the tunneling regime. The grain performs small oscillations 
around an average position, located between the the open region border and the lead. In the figures above, the open region is 
bounded by the vertical solid lines in the center of each junction. 
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FIG. 6. Plot of the position of the grain as a function of time for the bias voltage V = l.lVb and the parameters a = 6.4* 10~* 
and = 4.27 * 10~^. For this voltage the behavior of the system is on average very much like a static symmetric double 
junction. The jaggedness of the curve comes from the very different velocities of charged and uncharged grains. 
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1: The grain is charged from the near lead. 




2: Being positively charged, the grain is pushed 
towards the other lead, where two charges tun- 
nel off the grain. 




3: Being negatively charged, the grain is 
pushed back to the first lead, where the pro- 
cess starts over. 

FIG. 7. Schematic illustration of the charge transfer mechanism in the shuttle regime. The grain performs oscillations back 
and forth between the leads, loading and unloading two charges at each turning point. 
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Plot of the position of the grain as a function of time for the bias voltage V 

-3 



2.0Vb and the parameters a = 6.4*10 



FIG. 

and u)^ — 4.27* 10"'' . For this voltage an uncharged grain is everywhere in tunneling contact with both leads so that the charge 
transfer cycle illustrated in Fig. ^ is possible. 
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FIG. 9. Schematic illustration of the three parts of the average half-period for a shuttle cycle discussed in the text. Two 
different kinds of period times are illustrated. The further the grain moves towards the lead, the shorter the time T2 can be 
expected to be. 
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FIG. 10. Comparison between the current obtained by Monte Carlo simulations of the system shown in Fig. ^ (solid line) 
and the current as calculated by using the analytical expression in Eq. ( |2l| ) (circles). The parameters used in the simulation 
are: a = 6.4 * lO""* and lo"^ = 4.27 * 10"^. 
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